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Abstract: We conduct the first numerical simulations of lattice theories with exact super- 
symmetry arising from the orbifold constructions of ||l], ||, ^. We consider the Q = 4 theory 
in = 0, 2 dimensions and the Q = 16 theory in D = 0,2,4 dimensions. We show that 
the U{N) theories do not possess vacua which are stable non-perturbatively, but that this 
problem can be circumvented after truncation to SU{N). We measure the distribution of 
scalar field eigenvalues, the spectrum of the fermion operator and the phase of the Pfaffian 
arising after integration over the fermions. We monitor supersymmetry breaking effects by 
measuring a simple Ward identity. Our results indicate that simulations of AA = 4 super 
Yang-Mills may be achievable in the near future. 
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1. Introduction 

The study of supersymmetric theories on lattices has a long history - see the review Q and 
references therein. Recently there has been a resurgence of interest in the field with successful 
constructions of lattice theories which keep intact a subalgebra of the full supersymmetry 
algebra 

In this paper we will be concerned with specific discretizations of the Q = 4 and Q = 16 
supercharge Yang-Mills theories in a variety of dimensions. The lattice actions we employ 
were first derived using orbifold/deconstruction techniques in |l], |2|, ^. Recently, it was 
shown how to recover them by discretization of a twisted version of the target super Yang- 
Mills theories |^]^. Other proposals for lattice actions based on twisting can be found in 
El' li' l^' II' Ell- ^il^ language of the twisted constructions in this paper. 



^The connection between topological twisting and orbifold constructions had been anticipated earlier in 
nol hll] and has since been generalized in 
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The key feature of these actions is that, in addition to gauge invariance, they retain one 
or more exact super symmetries at non vanishing lattice spacing. Since this feature is already 
sufficient to pair each bosonic state with a fermionic state of the same energy it is expected 
that these models may flow to the target continuum theory with a minimum of fine tuning 
as the lattice spacing is sent to zero. 

In this paper, we show that the vacua of these lattice theories with U (N) gauge symmetry 
exhibit a non-perturbative instability associated to a runaway trace mode of the scalar fields. 
We show that this problem can be evaded if the gauge group is truncated to SU (N) at the 
price of a mild breaking of the exact super symmetry. In two dimensions and for the Q = 16 
theory we show that, nevertheless, supersymmetry is restored without fine tuning in the 
continuum limit. The situation in the two-dimensional Q = 4 theory is more problematic 
on account of the observed strong fluctuations in the phase of the Pfaffian that results from 
the fermion integration. These fluctuations lead to large statistical errors in all observables 
making it difficult to draw definite conclusions. 

In contrast we observe rather small phase fluctuations for the Q = 16 theory with gauge 
group SU{2) in both two and four dimensions for the small lattices used in this study. 

The organization of the paper is as follows; we first summarize the lattice actions with 
which we are concerned, describe some of the details of the simulation algorithms that are 
employed and then summarize our numerical results for the case of Q = 4 and Q = 16 
supercharge theories in dimensions zero, two and four. We end with a discussion and outlook 
for the future. 

2. Lattice actions 

2.1 Q = 4 theory in two dimensions 

The field content of the lattice theory comprises a multiplet of p-form fermions {tj, tp^,Xi2) 
distributed over sites and links of the lattice together with a complexified Wilson gauge link 
Ufj^. These fields transform under a scalar supersymmetry as follows: 



Qv 
Qd 




II. 











d 








Furthermore, as was shown in ||8| the action of the theory can be written in a Q-exact form 
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where k is a dimensionless bare coupling. The decomposition of the fermions into p-forms and 
the appearance of a scalar supersymmetry arises as a consequence of the twisting procedure 
which was described in detail in []l9| , [20t| . 

Notice that this supersymmetry is nilpotent making the supersymmetric invariance of 
the the lattice action in eqn. p.2| manifest. In this formulation the gauge links are non-unitary 
matrices of the form U^{x) = e"^M(^)+«^M(^) with ^ = 1,2. The imaginary parts of the gauge 
fields generate the scalar fields of theory in the continuum limit as was shown in Q ^. The 
gauge covariant difference operators appearing in eqn. |2.2| are defined by 



p{,+)/,(x) = ^/^(x)/,(x + e^) - /.(x)^/^(x + e,) (2.3) 

P^"Vm(x) = /m(x)Z7m W - - e^)Ui^ - e,) (2.4) 

where the unit lattice vectors are ei = (1,0), 62 = (0, 1) and the lattice field strength J^^i, is 
given by 

.F^, = Vl+^U,{x) = U^{^)U^{^ + /x) - U,{^)U^{^ + u) (2.5) 

In general the lattice fields are associated to links of the lattice and they transform corre- 
spondingly under gauge transformations: 

G(x)r?(x)Gt(x) 
G(x)V^^(x)Gt(x + e^) 

G(x)^/^(x)Gt(x + e^) 

G(x + e^)Z7^(x)Gt(x) (2.6) 



r/(x 



Notice that this choice of link and orientation for the twisted lattice fields maps exactly into 
their r-charge assignments in the orbifolding approach Indeed, the lattice actions derived 
bv the two methods are identical. 

The final lattice action takes the form 

5 = Tr (j^l^F^u + \ (vl~\) ' - X^uVl^^^u] - Vm) (2-7) 

where we have integrated out the auxiliary field d. 

It is also possible to consider dimensional reductions of this model to lower dimensions. 
For example we will show results obtained from simulations of the D = matrix model that 
results from this lattice theory after dropping all dependence on lattice coordinates x. The 
dimensionless lattice coupling k in reduced dimension D is given by 

where X = g^N is the 't Hooft coupling, (3 is the physical extent of the system and L the 
lattice length. 

^We use an antihermitian basis for the generators 
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2.2 Q = 16 theory in four dimensions 

Remarkably, the Q = 16 supercharge theory in four dimensions also possesses a Q-exact term 
in its action of precisely the same form as its two dimensional cousin: 



S = kQ 



(2.9) 



where the indices now run a,b = 1 ... 5 and the sixteen fermions of AA = 4 super Yang-Mills 
are now built from one scalar, five vectors and the ten components of the antisymmetric tensor 
Xab- The real parts of f/^, /i = 1 . . . 4 yield the usual four dimensional gauge field while 
and the remaining imaginary components of C/^ yield the expected six scalar fields of = 4 
super Yang- Mills |^]. 

The action of the nilpotent, scalar supersymmetry is the same as before 

QUa = 

QV'a = 
Q Xab = 

Qv = d 

Qd = (2.10) 
In addition a new Q-closed term is needed of the form 

5'closed = X] '^^ (^abcdeXdei^ + Sa + Bb + BcjV^c''^ + Gc) (2-11) 

X 

which is zero by virtue of an exact Bianchi identity satisfied by the lattice field strength 

eabcdeDi+^Tde = (2.12) 

Gauge invariance of this term requires = 0. This condition is satisfied if the basis 

vectors are taken to be ejj = 6^ with a = 1 ... 4 and 65 = (—1, —1, —1, —1) corresponding to 
their r-charge assignments in the orbifold construction j^. 

Again, the link assignments of the fields can be summarized by specifying the variation 
of the fields under gauge transformations 

??(x)^G(x)r?(x)Gt(x) 

V'a(x) ^ G(x)V'a(x)G(x + ea) 

Xafe(x) G(x-hea-Feb)xafe(x)G"^(x) 
Ua{^) ^ G(x)^/^(x)Gt(x + e,) 

Ua{x) ^ G(x + e,)Z7^(x)Gt(x) (2.13) 

(2.14) 
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The action of the covariant difference operators appearing in the action take the same form 
as for the Q = 4 theory eqn. 2.4. In the naive continuum hmit it is possible to show that 
this lattice theory reduces to the Marcus topological twist of AA = 4 super Yang- Mills (in 
modern parlance the GL- twist [^]) and hence in flat space is fully equivalent to the usual 
continuum theory. 

Again, the model can be reduced to lower dimensions by simply taking the fields to be 
independent of certain lattice coordinates and we will show data for the model in both zero, 



two and four dimensions. The bare coupling k is again given by eqn. 2. 



3. Simulation details 



Since both the Q = 4 and Q = 16 supercharge theories have such a similar structure they can 
be simulated using a single Monte Carlo code in which the lattice geometry, gauge group and 
number of supercharges are input as parameters. The dynamical twisted fermions are handled 
exactly and efficiently using the RHMC algorithm |2^. For completeness we summarize the 
main features of this algorithm as applied to the simulation of supersymmetric lattices here. 

If we denote the set of twisted fermions by the field ^' = (r/, ijj^, Xfj.i>) we first introduce a 
parallel pseudofermion field $ with action 



5PF = $"f(MtM)-3$ (3.1) 
where M = M(pl,'U'^) is the antisymmetric twisted lattice fermion operator given, for example. 



in eqn. 2.7 



P 



Integrating over the fields $ will then yield (up to a possible phase) the Pfaffian of the 
operator M{l/(,U^) as required. The fractional power is approximated by the partial fraction 
expansion 

, P 

where the coefficients {ai,f3i} are evaluated offline using the Remez algorithm to minimise 
the error in some interval {e,A). Typically we have used P = 15 which yields a fractional 
error of 0.00001 for the interval 0.0000001 1000.0 which conservatively covers the range we 
are interested in. 

Following the standard procedure we introduce momenta {puiPf) conjugate to the coordi- 
nates (U,^) and evolve the coupled system using a discrete time leapfrog algorithm according 
to the classical Hamiltonian H = Sb + Spp + puPu + One step of the discrete time 



^The antisymmetry is guaranteed if the fermion action is rewritten as the sum of the original terms plus 
their lattice transposes 



- 5 - 



update is given by 





5t 


opu = 




ntl^ — 
up(p — 




-h 


5U = 




5$ = 


5tp^ 






= 




6pq, = 







I\U 



where the forces fu and are given by 



fu 
h 



6S_ 

"SU 
5S_ 

"5^ 



(3.3) 
(3.4) 

(3.5) 
(3.6) 
(3.7) 

(3.8) 

(3.9) 
(3.10) 



and the bar denotes complex conjugation. Using the partial fraction expansion given in 



eqn. 3.2 these forces take the form 

p 

fu = ^oi. 

h 



i=l 



U 



5M 



6U 

p 

ao^ - ^ aiSi 

i=l 



6M 
olA 



where 



(MtM + Pi)si = $ 

ti = Msi 



(3.11) 

(3.12) 
(3.13) 



(3.14) 
(3.15) 



The latter set of sparse linear equations is solved using a multimass CG-solver |2J] which 
allows for the simultaneous solution of all P systems in a single CG solve. 

At the end of one such classical trajectory the final configuration is subjected to a stan- 
dard Metropolis test based on the Hamiltonian H. The symplectic and reversible nature of 
the discrete time update is then sufficient to allow for detailed balance to be satisfied and 
hence expectation values are independent of 5t. After each such trajectory the momenta are 
refreshed from the appropriate gaussian distribution as determined by H which renders the 
simulation ergodic. 

In the work reported here we have employed only periodic boundary conditions which 
preserve supersymmetry. Thermal boundary conditions are also interesting as they allow for 
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exploration of dualities between string and gauge theory 26, 2^, 2£, 3C] and dynamical 



supersymmetry breaking |31, 32 1. 



Our measurements concentrate, in part, on local observables such as the eigenvalues of 
l4ji{x)L{fj,{x) and the bosonic action < Sb{1^) >• The former yields, in the continuum limit, the 
distribution of eigenvalues of the scalar fields Z/Z^Z/Z^j—/ = 2i?^+. . . and hence gives information 
on the quantum moduli space. The latter observable is related to a exact supersymmetric 
Ward identity and can be evaluated analytically which provides both a useful check on our 
code and measures the magnitude of supersymmetry breaking effects. The analytic argument 
that determines the value of < Sb > is simple. Consider first the Q = 4 supercharge theory 
and write down an expression for the mean action < 5 > 

< 5 >= -^1^ =< QA >= (3.16) 

OK 

where the last result follows from Q-exact nature of the twisted action and shows that the 
vanishing mean action is the consequence of a simple Q-Ward identity. This argument needs 
a minor modification for the Q = 16 supercharge theory which contains also a Q-closed term. 
However, it is straightforward to show that a simple rescaling of the field X/xi/ — *■ y/^Xfiu 
removes the k dependence of this term so that once again = as a consequence of 

Q-sup er symmetry^ . 

The result < S >=< Sb > + < Sp >= can be translated into an exact result for the 
bosonic action since the fermions appear quadratically in the action and hence < Sp > can 
be evaluated by a simple scaling argument. In the case of Q = 4 supercharges one finds 

k<Sb>= ^NgV (3.17) 

where Nq is the number of generators of the group and V is the number of lattice points. 
One might have naively expected a factor of 4 representing the four twisted fermions of the 
Q = 4 theory rather than the factor of 3 that is present in this expression - the discrepancy 
arises as a consequence of integrating out the auxiliary field d which effectively removes the 
contribution of one fermion to < Sb >■ 

The Q = 16 theory is a little more involved. To evaluate < Sp > by a scaling argument 
requires an additional rescaling of both t] and ip^ by a factor of ^/K. This results in an 
additional multiplicative factor of the measure. The final result for < Sb > is 

then 

k<Sb>= ^NgV (3.18) 

where the factor arising in the numerator is now composed from 9 = 16 — 6 — 1. 

A faster way to derive these results relies simply on the coupling constant independence 
of the free energy; the bosonic action can then be evaluated in the weak coupling limit where 



^Another way to see this is to realize that the partition function for periodic boundary conditions is just 
the Witten index and hence does not depend on the coupling constant 
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the theory is quadratic in the bosons and equipartition holds; the bosonic action then simply 
counts the number of degrees of freedom. 

For the small systems we have examined in this paper we have also measured the Pfaffian 
and the spectrum of the fermion operator M(L{,U). The Pfaffian computation is carried 
out by using a variant of Gaussian elimination with full pivoting to transform the 2n x 2n 
dimensional antisymmetric matrix M into the canonical form 








Ai 





. 


..\ 




-Ai 








. 















A2 . 












-A2 


. 




v 











. 





(3.19) 



Then 



Pi{M) = llXi 



1=1 



(3.20) 



As can be seen by examining eqn. |3.1| our simulations generate the phase quenched ensemble 
defined by \Pf{M)\. As usual we can always compensate for neglecting any phase e*"^^^ by 
re- weighting all observables 0{U) by the phase factor according to the simple rule 



<0>-- 



< 0{U)e 



ia{U) 



>a=0 



Where necessary we have reweighted our results accordingly. 



(3.21) 



4. Zero dimensions 

4.1 Vacuum instability for U{N) theories 

The construction we have described is strictly valid only for U{N) theories. This can be 



seen in a variety of ways; taking the trace of the first line of the Q- variations in eqn. 2.1C 
is inconsistent if the fields are restricted to the traceless generators. Similarly the difference 
operators employed in the lattice action when applied to a traceless field generically yield 
a field with non-zero trace. Thus we initially focus on simulations of the U{N) theories. 
Unfortunately we will see that these theories have a non-perturbative instability; the trace 
mode of the scalar fields rapidly runs off to (negative) infinity. 

Figure ^ shows this explicitly with a plot of the Monte Carlo history of the eigenvalues 
oi ul,{x)Uf_i{x) for an SU{2) model with Q = 4 supersymmetries reduced to zero dimensions. 
Clearly one of the eigenvalues is driven to zero after a finite simulation time. Using the 
representation 

U^{x) = c^/^W+^^mW (4.1) 
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Q=4 D=0 U(2) m=0.1 




Figure 1: Eigenvalues of U^Uf^ versus MC time. Gauge group SU{2), Q — A and D = 



m 


Sb 


Qcxact 




0.01 


2.65(6) 


6.0 


0.45(2) 


0.1 


2.40(8) 


6.0 


0.57(6) 


0.5 


2.99(7) 


6.0 


0.38(2) 



Table 1: Mass dependence of ^TTUj^Uf^ 



it is easy to see that det(M(a;)) = e^^^^°^^^ where B'^{x) is the trace component of the 
scalar field^. The presence of a zero eigenvalue leads to a vanishing determinant and clearly is 
associated with a limit in which the trace mode B^. — > — oo. In this situation the fluctuations in 
the gauge links are never small and no limit exists in which the lattice model approximates the 
correct continuum theory. Furthermore, we have observed that no simple scalar mass term as 
advocated in is able to save the situation. Table |l] shows the value of < -^Tr u}i{x)Uf^{x) > 
obtained from a series of runs in which the lattice action is supplemented with the additional 
term 

' ^ ' (4.2) 



The first column gives the mean bosonic action as compared, in the second column, to 
its exact value as predicted by eqns. 3.17 and 3.18i Clearly, the fields deviate strongly 
from the unit matrix for all values of the mass parameter m and furthermore supersymmetry 
is also badly broken. In all these cases we also observe a zero eigenvalue in M independent of 
the mass. These conclusions survive for other gauge groups U{N), values of the supercharge 
Q and dimension of the model. 



Clearly one zero eigenvalue is the minimum required for consistency with a vanishing determinant - we 
have also encountered runs where some or all of the UjjA^ matrices have more than one zero eigenvalue 
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It is straightforward to see why this occurs; the bosonic action evaluated on an arbitrary 
{l{^{x)} configuration is generically lowered by a global shift in the trace mode -B^(x) 
B^{x) — c with positive c^. Explicitly, 

SBil(^{x)e-'') = e-^'SB(p(^{x)) (4.3) 

since e""^^ commutes with all terms in the bosonic action (each of which contains a product 
of four matrices). Unlike in the continuum the penalty from the mass term is finite even 
in the limit —oo and entropic effects appear to drive the system far from the point 

S° = 0. Notice that this effect occurs only for the exponential parametrization of the complex 
matrices used here and given in eqn. |4.1| . The prescription utilized in the original orbifold 
constructions used instead the decomposition 

U^ = I + A^ + iB^ (4.4) 

which is not subject to this same instability. 

One way to evade these problems immediately presents itself; simply set the trace mode 
to zero by truncating the lattice model to the special unitary group SU{N). Notice that 
this truncation preserves supersymmetry if the (complexified) gauge fields lie infinitessimally 
close to the identity hl^{x) = I + aA^{x) + . . . which should happen in the continuum limit. 
Presumably the restriction to SU{N) is also irrelevant in the large N limit. However, for 
non-zero lattice spacing and finite N it will lead to supersymmetry breaking effects which we 
examine in the next section. 

4.2 Vacuum structure for SU{N) theories 
4.2.1 Q = 4 supercharges 

Once the truncation to SU (N) is carried out we observe no instability; the eigenvalues oiUliU^ 
cluster around unity for any gauge group, dimension or number of supercharges. Further- 
more the width of this eigenvalue distribution decreases as the bare lattice coupling increases 
which is a necessary condition for the lattice model to approximate the target theory in the 
continuum limit. 

Nevertheless, after this truncation we might expect to see supersymmetry breaking effects 
and we have probed for this by examining the bosonic action. As we have argued earlier the 
expectation value of this term can be derived exactly as a consequence of the twisted scalar 
supersymmetry. Deviations from the exact value hence measure supersymmetry breaking. 
Table ^ shows results from a simulation of the Q = 4 model with SU{2) gauge group in 
the matrix model limit^. While a small, statistically significant deviation is seen for strong 
coupling this disappears for larger k. The limit k — > oo corresponds to the continuum limit 
a — >■ for theories in 1 < D < 4. A similar pattern seen in table ^ for the same model with 
group SU{3). 

®For a classical vacuum state corresponding to constant diagonal matrices the bosonic action remains zero 
under this shift 

'^The Pfaffian is real positive definite for all A*' in this limit and hence no issue of reweighting occurs 



- 10 - 





<Sb> 


Qcxact 


1.0 


4.40(2) 


4.5 


10.0 


4.47(2) 


4.5 


100.0 


4.483(15) 


4.5 



Table 2: Bosonic action for SU{2) Q = 4 model at several couplings 



K 


< Sb> 


Qoxact 


1.0 


11.71(2) 


12.0 


10.0 


11.98(3) 


12.0 


100.0 


11.98(4) 


12.0 



Table 3: Bosonic action for SU (3) Q = 4 model at several couplings 



Q=4 D=0 SU(2) I 

2.5 - 




5 10 15 20 25 30 

lambda 

Figure 2: Probability distribution of eigenvalues of UjjJfj, — I for Q = 4 and SU (2) and D = 

Having verified that the breaking of supersymmetry is indeed small after the truncation 
from U (N) to SU {N) we now turn to the scalar field eigenvalues. Figure |2| shows a plot of 
the probability distribution of the eigenvalues of UjiU^ — I for SU (2) gauge group. Classically 
the model contains flat directions corresponding to constant diagonal U matrices. One might 
worry that the presence of these flat directions might render the path integral defining the 
quantum theory ill defined but these numerical results, like earlier matrix model studies 

|3^ , confirm that the scalar fields are localised around the origin in moduli space and 
the partition function is finite*. The asymmetry in the distribution is a cut-off effect; the 

*As has been seen before the distribution of scalar field eigenvalues is observed to possess A'^ peaks for 
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Figure 3: Probability distribution of eigenvalues of the fermion operator for Q = 4, SU{2) and D = 



lattice theory, unlike the continuum theory, is not invariant under a change of sign of the 
scalar fields. Indeed, the entire region — oo < < is mapped onto the finite segment 
— 1 < UjiUf^ — 1 < 0. The asymmetry in the distribution can be hence be used to quantity 
the magnitude of lattice artifacts. Since the dimensionless eigenvalues contain a factor of the 
lattice spacing the tails of the distribution contract and are less sensitive to the hard cut-off 
A = — 1 as the lattice spacing is reduced. 

More importantly, notice that the effect of the classical flat directions is still visible in the 
presence of a power law tail in the distribution extending out to large (positive) eigenvalue. 
Both theoretical arguments and numerical results suggest the power p = 3 independent 



of for this model [34|. Such a power behavior would yield a logarithmically divergent 
value for < > which means that in any Monte Carlo simulation the scalars spend an 
appreciable amount of time far the origin in field space. In the background of one of these 
field configurations the bosonic action will necessarily develop a near zero mode corresponding 
to translations along this flat direction. By supersymmetry we expect the fermion operator 
must also then develop a near zero mode. This can be seen explicitly in figure |3| which 
shows the distribution of the absolute value of the fermion eigenvalue for this system. An 
enhancement at small eigenvalue is seen consistent with the previous argument. 

4.2.2 Q = 16 supercharges 

In tables g | we examine the bosonic action for the Q = 16 supercharge model in the matrix 
model limit. Notice that in the case of Q = 16 supercharges the fractional deviation of 
the expectation value from its exact supersymmetric value is smaller for SU (3) than SU{2) 
at fixed lattice coupling which is consistent with the breaking effect vanishing in the limit 
N ^ oo. The same effect was not visible however for Q = 4 supercharges. 

group SU{N) 
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<Sb> 


cexact 


1.0 


13.67(4) 


13.5 


10.0 


13.52 


13.5 


100.0 


13.48(2) 


13.5 



Table 4: Bosonic action for SU{2) Q = 16, D = model at several couplings 





< Sb> 


Qcxact 


< cos a > 


1.0 


36.01(20) 


36.0 


0.67(5) 


10.0 


35.95(18) 


36.0 


0.47(12) 


100.0 


35.92(10) 


36.0 


0.47(13) 



Table 5: Bosonic action for SU(3), Q = 16, D = model at several couplings 



10 



6 



4 



2 
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, 1 1 1 1 1 1 1 , 1 

Q=16D=0SU(2) - 




i/^*^^^ 1,1,1, 



Figure 4: Probability distribution of eigenvalues of UjjUf^ — I ioi Q ~ 16, SU{2) and D ~ 

Additionally notice that table ^ contains a fourth column corresponding to the average 
value of the cosine of the Pfaffian phase (the average sine being always consistent with zero) . 
While the Pfaffian for the Q = 16 supercharge model is generically complex it is possible 
to show that it is real positive definite for the SU{2) model in the matrix model limit. In 
the case of SU{3) the Pfaffian in zero dimensions is real but not necessarily positive definite 
and hence the bosonic action is reweighted with the sign of the Pfaffian, as discussed in the 
previous section. 

The probability distribution for ujiU — I for Q = 16 and SU{2) is shown in figure ^. 
While it is qualitatively similar to the Q = 4 case it should be clear that the localization 
around the origin is more dramatic than for Q = 4 ~ indeed, theory and simulation point to 
a larger value of the power law exponent p = 15 governing the tail of the distribution. Since 
now the scalar fields do not penetrate far down the classical flat directions one would expect 
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1.5 

-b 

B 1 
0.5 



"0 1 2 3 4 5 6 7 

I lambda I 

Figure 5: Probability distribution of eigenvalues of the fermion operator for Q — 16, SU{2) and 
£1 = 



K 


< Sb >" 


<Sb> 


cexact 


< cos a > 


8.0 


17.15(1) 


17.56(3) 


18.0 


-0.24(1) 


18.0 


39.23(2) 


41(6) 


40.5 


-0.06(2) 


32.0 


70.61(4) 


65(5) 


72.0 


-0.014(6) 



Table 6: Observables for SU{2) Q = 4 model in £> = 2 



to see correspondingly fewer near zero modes in the fluctuation spectra of either the bosonic 
or, by supersymmetry, fermionic operator in this case. Figure |5| shows that indeed this is 
the case - a gap appears to open up in the probability distribution of the magnitude of the 
fermion eigenvalue close to the origin. 

We will see that these qualitative features survive in the two dimensional case to which 
we now turn. 

5. Two dimensions 
5.1 Q = 4 supercharges 

We have also begun an investigation of the Q = 4 model with gauge group SU{2) in two 
dimensions. Again, we have focused on the bosonic action, the distribution of scalar eigen- 
values and the spectrum of the fermion operator. In these simulations we hold fixed the 
continuum dimensionless 't Hooft coupling = 0.5 as defined in eqn. which leads to 
a lattice coupling that grows like the square of the lattice length L. Table |^ shows results 
from simulations on lattices of size L = 2,3,4 for gauge group SU{2). The second column 
corresponds to the phase quenched approximation in which the data is not reweighted by the 
Pfaffian phase. The third column corresponds to the reweighted action. Notice that the phase 
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fermion eig Q=4 D=2 SU(2) 




2 3 4 5 6 7 

I lambdal 



Figure 6: Probability distribution of absolute fermion eigenvalue for Q = 4:D = 2L = 2 and SU{2) 



quenched numbers are relatively close to their exact values. However, while there is a hint 
from the data that reweighting pushes them closer to those exact values, it should be clear 
that the statistical errors increase rapidly with lattice size rendering it impossible to make 
meaningful measurements on the larger lattices. This is confirmed by looking at the average 
of the cosine of the Pfaffian phase shown in the fifth column. It is statistically consistent with 
zero on the larger lattices. 

Actually, the observation that the phase factor vanishes in this theory is quite interesting; 
since we employ periodic boundary conditions the measured expectation value yields the 
Witten index for the theory: 



— < e ^ phase quenched 



(5.1) 



Thus our data is consistent with a vanishing Witten index - a necessary condition for su- 
persymmetry breaking. Furthermore, as in the matrix model case, we have observed a large 
number of rather small fermion eigenvalues - figure ^ shows the distribution of the absolute 
value of the fermion eigenvalue for the SU{2) model on a lattice with L = 2. The existence of 
many small magnitude eigenvalues is highlighted by the scatter plot of figure ^ showing the 
real and imaginary parts of all fermion eigenvalues with magnitude |A| < 1.5 obtained from 
a sample of 1000 configurations for L = 2. Notice that the eigenvalues, while concentrated 
in a band along the imaginary axis, nevertheless have a non-zero density, by virtue of the 
Yukawa couplings, over the entire plane including the region around the origin. As we will 
see later this last feature is not seen for Q = 16 supercharges where we observe essentially no 
eigenvalues in the vicinity of the origin. 

For Q = 4 it is possible that these near massless states could play a role as Goldstino 
modes associated with dynamical supersymmetry breaking which has been conjectured to 



happen for the Q = 4 theory in low dimensions [36|. 

The scalar eigenvalue distribution is similar to its cousin in zero dimensions and is shown 
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Figure 7: Eigenvalues of the fermion operator for Q = 4, SU{2) and D ^ 2 
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Figure 8: Probability distribution of eigenvalues of UjjJf^ — / for Q = 4, SU{2) and D = 2 



in figure g and shows again that the scalar distribution possesses a tail to large eigenvalue. 
Furthermore, we have observed that there is a correlation between the occurrence of small 
fermion eigenvalues and the presence of large scalar fields associated with a significant com- 
ponent of the scalar fields along the classical flat directions. 

Our overall conclusion is that practical simulations of the Q = 4 supercharge theory in two 
dimensions with supersymmetry preserving periodic boundary conditions will be extremely 
hard due the rapid fluctuations in the phase of Pfaffian - a problem which was first highlighted 
in Furthermore, it is possible that these phase variations are associated with a dynamical 
breaking of supersymmetry as argued for in |36]. Our numerical results support a vanishing 
Witten index and the presence of massless fermions - necessary conditions for supersymmetry 
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ccxact 


< cos a > 


8.0 


53.26(6) 


53.26(6) 


54.0 


0.999997(1) 


18.0 


120.1(2) 


120.1(2) 


121.5 


0.999995(1) 


32.0 


214.7(4) 


214.6(3) 


216.0 


0.999994(3) 



Table 7: Observables for SU{2) Q = 16 model in D = 2 




Figure 9: Probability distribution of ferniion eigenvalues for Q = 16, D = 2, L = 2 and SU{2) 



breaking. If so these results contradict the numerical results presented in |31]. A spontaneous 
breaking of supersymmetry would also invalidate the theoretical arguments given by Matsuura 
[^] showing that the orbifold theories have zero vacuum energy. However the latter result is 
derived using semiclassical exactness which is invalid if the Q-supersymmetry is spontaneously 
broken by non-perturbative effects. Clearly, further work is required to resolve this question 
unambiguously. 



5.2 Q = 16 supercharges 

The results for the bosonic action for the Q = 16 model are shown in table |^ and correspond 
to simulations on lattices of size L = 2,3,4 and gauge group SU{2). Perhaps the most 
striking result is that the Pfaffian phase is very close to zero and remains so as the lattice size 
is increased and the continuum limit approached. In this case the reweighting procedure is 
irrelevant as can be seen by comparing the phase quenched and full expectation values shown 
in table ^. 

We believe that this suppression in the Pfaffian phase is related to the phenomena we 
observed in zero dimensions; the scarcity of small fermion eigenvalues. A plot of the distri- 
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Figure 10: Eigenvalues of the fermion operator for Q — 16, SU{2) and D — 2 
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Qcxact 


< cos a > 
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211.2(2) 


211.2(2) 


216.0 


0.9945(3) 


3 


1072.8(10) 


1075.0(35) 


1093.5 


0.955(6) 



Table 8: Observables for SU{2) Q = 16 model in D = 4 at A = 0.5 



bution of the absolute fermion eigenvalue for this theory is shown in figure ^ and indeed a 
gap appears to open up in the spectrum. This conclusion is reinforced when we examine a 
scatter plot of the fermion eigenvalues in the complex plane given in figure |lO[ Remarkably, 
and in contrast to the Q = 4 model, the eigenvalues are excluded from a region around the 
origin and concentrate along the imaginary axis. It is reasonable to conjecture that it is the 
small eigenvalues that control the phase; after all in the limit in which the eigenvalues are 
confined to the imaginary axis an eigenvalue must flow from positive to negative values for 
a sign change to occur. Such a transition would require a zero eigenvalue to arise which we 
observe to be highly unlikely. Thus phase fluctuations are suppressed. 

Furthermore, the presence of a gap in the fermion spectrum correlates to a rather rapidly 
damped scalar eigenvalue distribution for the Q = 16 model as we saw in zero dimensions. 



This distribution is shown in figure 11 



6. Four dimensions 

Finally we report on our preliminary simulations of the Q = 16 supercharge theory in four 
dimensions. Table ^ shows the bosonic action and cosine of the Pfaffian phase for lattices 
with size L = 2, 3 at fixed 't Hooft coupling A = 0.5 (the data corresponds to 6000 and 1000 
configurations for L = 2 and L = 3 respectively) 
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Figure 11: Probability distribution of UjjJfj. — I eigenvalues for Q — 16, D — 2, L — 2 and SU{2) 
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ccxact 
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212.5(3) 


211.5(5) 


216.0 


0.9964(2) 


3 


1081.9(17) 


1080.5(45) 


1093.5 


0.983(2) 



Table 9: Observables for SU{2) Q = 16 model in D = 4 at A = 0.25 



While supersymmetry breaking effects are visible at 0(1)% the Pfaffian phase is small 
and reweighting offers a reliable way to deal with the fluctuations at least for small lattices^. 
It remains to be seen whether this is true as the lattice size is increased but we see these 
results as encouraging (after all we have lost control of the phase already with small lattices 
for the Q = 4 model). 

For comparison, table |^ shows the same quantities for 't Hooft coupling A = 0.25. Notice 
that as we approach weak coupling and smaller lattice spacings the bosonic action moves 
towards its exact super symmetric value as expected. Furthermore, the phase fluctuations 
also decrease. The question of whether reweighting will be practical on larger lattices will 
hence depend on which trajectory Abare(-^) we must follow in the (A, L) plane to approach 
the continuum limit. 

The scalar and fermion eigenvalue distributions are shown in flgures 12 and ^ and look 
qualitatively similar to their two dimensional cousins. 

Again, the scatter plot of fermion eigenvalues also indicates that eigenvalues are repelled 
from a region around the origin as for the Q = 16 theory in two dimensions. 



®Tlie larger error in tiie reweighted numbers for L = 3 reflects only the decreased statistics available in tiiat 
case wliere tiie Pfaffian is only computed every tentii measurement 
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Figure 13: Probability distribution of scalar eigenvalues ior Q — 16, D = 4:, L — 2 and SU{2) 



7. Discussion 

In this paper we have presented numerical results from simulations of a variety of Yang- 
Mills theories with extended super symmetry. The lattice actions that are employed can be 
derived either with orbifold methods [|l|, ^, |^ or via geometrical discretization of a twisted 
version of the target theory Remarkably, they possess both gauge invariance and one 
more exact supersymmetries at non-zero lattice spacing. These lattice theories are important 
both at a conceptual and practical level; they offer the possibility of a rigorous definition of 
the continuum gauge theory, and through numerical simulation may offer up new ways to 
extract non-perturbative information on that gauge theory. In this paper we have focused on 
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Figure 14: Eigenvalues of the fermion operator for Q — 16, D = i and SU{2) 

the latter question; specifically, are these lattice theories amenable to Monte Carlo simulation 
using the tools and techniques of lattice gauge theory ? 

We have shown that the U{N) theories generically suffer from a vacuum stability problem 
- the theory is defined in terms of complexified Wilson gauge links and these develop one 
or more zero eigenvalues under quantum corrections. This instability is associated with the 
trace mode of the would be scalar fields in the lattice theory. As a result the link fields are 
driven a long way from the identity for any value of the lattice spacing rendering invalid the 
correspondence between the lattice model and the target continuum theory. 

However we have shown that this instability of the lattice action is avoided if the theory 
is truncated to SU{N)^^. Furthermore, by measurement of a simple super symmetric Ward 
identity we have presented evidence that the associated breaking of supersymmetry is small 
and decreases as the continuum limit is taken. 

A second potential problem arises however; integration over the fermions in these models 
generically leads to a complex Pfaffian (which can be reduced to a determinant for Q = 4 
models). Our simulations are necessarily performed in the phase quenched ensemble where 
this phase is ignored. In principle, expectation values can be computed in the full ensemble 
by a reweighting procedure. However, this is only practical if the phase fluctuations are small. 
In the case of the Q = 4 theory we find the average phase factor approaches zero rapidly as 
the continuum limit is taken and the statistical errors in the reweighted observables grow 
uncontrollably. This seems to rule out the possibility of using these actions to study the 
Q = 4 theory at least for zero temperature (periodic temporal boundary conditions for all 
fields). 

Intriguingly, the vanishing expectation value of the phase factor < e*"^ > phase quenched = 
may be given a physical interpretation - it corresponds naively to a vanishing Witten index 



'this truncation also has the merit of removing an exact zero mode of the fermion operator 



for the model. Of course the notion of a Witten index as a signed sum over classical vacua 
is somewhat delicate in theories with extended supersymmetry which possess a continuum of 
such vacua. Nevertheless, a vanishing Witten index is a necessary condition for supersymme- 
try breaking and our analysis of the fermion spectrum shows that indeed this theory contains 
one or more near zero modes which could play the role of a Goldstino associated with dynam- 
ical breaking of supersymmetry. An argument for such a breaking has been independently by 
Hori et al for the Q = 4 theory in two dimensions. 

Furthermore, we have observed that the appearance of small fermion eigenvalues is corre- 
lated with large excursions of the scalar fields along the classical directions. In the background 
of such a vacuum configuration the operator describing small fluctuations of the scalars de- 
velops a small eigenvalue corresponding to motion along the flat direction. This eigenvalue 
is mirrored in the fermionic sector because of supersymmetry. The question of whether this 
effect can lead to supersymmetry breaking is then tied to the question of how frequently, in 
the context of a Monte Carlo simulation, the scalars probe these flat directions which, in turn, 
is measured by the probability distribution of scalar field eigenvalues. For Q = 4 we observe 
that this distribution P{X) shows a slow power law decay which has been estimated to vary 
as for large A [33, 34]. This could supply the effect we have argued for; the variance of A 
then diverges logarithmically and the typical scalar field configuration wanders appreciably 
away from the origin and leads to a non-zero density of small fermion eigenvalues. 

At first sight the strong phase fiuctuations visible in the lattice Q = 4 theory are surprising 
since a proof exists that the continuum Pfaffian is real positive semi-definite and one might 
have therefore expected the phase fluctuations to be suppressed for small enough lattice 
spacing. This does not appear to be the case. However, we believe there is no inconsistency; 
if supersymmetry breaks spontaneously as we claim, the Pfaffian will be zero on the important 
continuum field configurations dominating the path integral, which is then consistent with 
the vanishing of the expectation value of the phase factor seen in the lattice simulations. 

In the case of the Q = 16 supercharge model the corresponding scalar distribution falls 
much more rapidly with eigenvalue. Correspondingly, we observe that the fermion spectrum 
shows a complete absence of near zero modes and the Pfaffian phase is typically small and can 
handled, at least for these small lattices and for SU{2), by reweighting. This is naively rather 
unexpected - afterall the Pfaffian is generically complex for the Q = 16 model. However, 
as we have argued previously, the average phase factor in the phase quenched ensemble is 
nothing more than the Witten index of the model which is a topological invariant. It can 
hence be evaluated exactly in the semi-classical limit along the lines described in provided 
supersymmetry remains unbroken which is thought to be the case for the Q = 16 model. Such 
a calculation shows that W = 1. Thus topological arguments would suggest that it should 
be possible to handle the Pfaffian phase for the Q = 16 theory using reweighting techniques 
even though the Pfaffian of the theory is generically complex. This is a remarkable result but 
quite consistent with the results of our simulations. 

The absence of small fermion eigenvalues together with the seeming relative unimportance 
of Pfaffian phase fiuctuations lead us to feel cautiously optimistic about the feasibility of 
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conducting large scale simulations of the Q = 16 supercharge theory in four dimensions in 
the near future. Such simulations are of course very interesting from the point of view of the 
AdSCFT correspondance and we hope to report on results from such simulations in the near 
future m. 
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